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1 Inroduction 

The n x n cyclic nonadiagonal matrix is as in the following form: 
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where n > 8. 

This type of matrix appears in many areas such as engineering 
applications. The determinants and the inversions of these matrices are 
usually required. Many algorithms are composed by using the LU 
factorization for the periodic tridiagonal, pentadiagonal and the cyclic 
pentadiagonal and heptadiagonal matrices pQ- [5] 

A new recursive symbolic algorithm for inverting general periodic 
tridiagonal and anti-tridiagonal matrices are studied in [2 . The authors 
compose a new symbolic algorithm for the inverses of the periodic 
pentadiagonal matrix and the periodic anti-pentadiagonal matrix is obtained by 
using it in [3]. With some restrictive conditions, algorithms for the 
inverses of the tridiagonal and pentadiagonal matrices are given in [5] . In [3] it is 
presented that a new computational algorithm to evaluate the determinant of 
the tridiagonal matrix with its cost. In 6 an expression of the 
characteristic polynomial and eigenvectors for pentadiagonal matrix is obtained 
and an algorithm to compute the determinant of the pentadiagonal matrix is 
presented. 

In this paper, we extend the work presented in [T]. In the second section we 
obtain the Doolittle LU factorization of the cyclic nonadiagonal matrix. Then 
by using the elements of the last six columns, the elements of remaining (n — 6) 
columns are found and the inverse matrix of the cyclic nonadiagonal matrix is 
obtained. In the last section a numerical example is given. 
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2 Main Result 



In this section, t which is only a symbolic name is chosen as a parameter. Then 
the determinant and the inverse of the cyclic nonadiagonal matrix K in (JlJ are 
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where 
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The determinant of the cyclic nonadiagonal matrix K is computed as in the 
following: 

n 

det(K) = Y[c, 

1=1 

When K is nonsingular, let 

K 1 — (Sij)i<i.j< n = (Ci, C2, . . . , C r , . . . , C n ) 
and Ci is the «th column of K~ x where 1 < i < n. Notice that 

for r = 1, 2, n. We can write C r as follows: 

C r = (Ci, C2, . . . , C r , . . . , C n )E r (2) 

where E r — (8i r , ^2r, • ■ • , $rr, ■ ■ ■ , 3in) T , r — 1,2, ...,n (6ij is the Kronecker 
symbol). By using ^ we obtain 

KQ = E t (3) 

for i = n, n — 1, n — 2, n — 3, n — 4, n — 5. 

Now, an algorithm can be composed by using the last six columns of K^ 1 . 
By using the LU factorization and ([3]), the components of the last six columns 
are computed as in the following: 
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S n -l,n-2 — -^-^{ — k n -2 — Vn-lS n ,n-2) 

S n —2,n—2 = c _ 2 (1 — ^>n-2 $71- 1,71—2 — ll 7i-2<SVi :n _2) 

S n ,n-3 — — ^n-3 + ^7i-2./n-2 + ~ ^n-l^n-2 /n-2) (4) 

Sn— l,n-3 — ^T"^ - ^"-3 + — Vn-lS n ,n-3) 

S n -2,n-3 — -^—;{ — fn-2 — W n -2iSn- 1,71-3 — ^ l n-2'S'n,n-3) 

SW— 3,n-3 — STjO- — e n-3>S7i-2,7i-3 — ^n-S'S'n-l.n-S — Wn-S^n.n-s) 
Sn,7i-4 = + hn-3fn-3 + ^n-2.9n-2 _ hn-lfn-lfn-3 

+ /l„_lfc„_4 — h n -\k n -?>fn-?> — h n -lkn-29n-2 + ^71-1^71-2/71-2/71-3) 

Sn-1,71-4 = ^( — ^71-4 + fcn-3/n-3 + k n -29n-2 — fcn-2/n-2/n-3 

Wji— l<5n,7i— 4) 
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S n -2,n-4 = 7~~;(~9n-2 + fn-2.fn-3 — W„_2>S'n-l,n-4 — V n-2<Sn,n-4) 
S n -3,n-4 = — — e„_35„_2,n-4 
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for j = n, n — 1, n — 2, n — 3, n — 4, n — 5 and « = n — 6, n — 7, 1 

S'iJ = (ej/Si+ij + PiSi + 2,j + TiSi+3j + ZiSi + 4.j + WiS n -l,j + ViS n .j) (10) 

The elements of the remaining (n — 6) columns can be calculated using the 
fact that K~ X K = I n where J„ is an n x n identity matrix. Then 
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C n -6 — {E n -2 — M n -5C n -5 — A n -4,C n -4, — a n -3C n -3 (H) 

Zn— 6 

— d n -2C n -2 — b n _iC n -i — B n C n ) 

C n -7 = (-En-3 — Af„_gC„_6 — A„_sC n _5 — fl„-4C n _4 (12) 

— rf„_3(7„_3 — b n - 2 C n -2 — B n _iC n -i — N n C n ) 



Cj = —{Ej+A-M j+ iCj + i-A j+ 2Cj + 2-aj + 3Cj + 3-dj + iCj + A (13) 
—bj +5 Cj + s — Bj +& Cj+6 — Nj + - l Cj + 'j — Rj + $Cj + %) 

where j = n — 8, n — 9, . . . , 1 and Zi ^ for i = 1, 2, . . . , n — 6. 

Thus, the algorithm for the n x n cyclic nonadiagonal matrix is given as 

Input: n is the order and di,di, Ai, Mi, Zi,bi, Bi, iVj, Ri are the entries of 
the cyclic nonadiagonal matrix. 

Output: Inverse matrix K^ 1 — (Sij)i<ij< n . 

Stepl:If Zi = for any i = 1, 2, . . . , n — 6 set Zi = t. 

Step2: If Ri = for any i = 6, 7, . . . , n set Ri = t. 

Step3: Set c\ = d\, if c\ = then ci = t, 
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c 2 = d 2 - f 2 ei, if c 2 = then c 2 



h = — (&3-ff3ei) 

54 = — (B 4 -a 4 e\) 

e 2 = a 2 - /2-P1 
P 2 = A 2 -/ 2 T! 
T 2 = M 2 -f 2Zl 

k 2 = (fcid) 
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w 2 = -/2W1 

/i 2 = — (A„-/nei) 
c 2 
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C3 = <fe - 03-Pi - /3e 2 , if c 3 = then c 3 = t 

h = —{h- olaP\ - g 4 e 2 ) 
C3 

e3 = 03 - 53^1 - /3-P2 
^3 = M - g 3 z x - f 3 T 2 
T 3 = M 3 -f 3 z 2 

k 3 = {kiPi + k 2 e 2 ) 

c 3 

w 3 = -g 3 wi - f 3 w 2 

h 3 = (hiPi + h 2 e 2 ) 

C3 

v 3 = -93V1 - hv 2 

c 4 = d 4 — ct 4 T\ — g 4 P 2 — f 4 e 3 if c 4 = then c 4 = t 

e 4 = a 4 - a 4 zx - g 4 T 2 - f 4 P 3 
Pi = Ai- g 4 z 2 - /4T3 
T 4 = M 4 -f 4 z 3 

k 4 = - — (fciTi + k 2 P 2 + k 3 e 3 ) 
c 4 

w 4 = —a 4 wi — g 4 w 2 — f 4 w 3 

h 4 = -—(h 1 T 1 + h 2 P2 + h 3 e 3 ) 
c 4 

v 4 = -a 4 vi - g 4 v 2 - f 4 v 3 



Step4: For i = 5, 6, . . . , n — 3 do 
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Step5: i = 5, 6 . . .n — 6 

fc» = -j:(h-4Zi-4 + ki- 3 Ti_ 3 + ki- 2 Pi-2 + fci_iej_i) 
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W„- 5 = - 7n-5 W n-9 - «n-5^n-8 - 9n-5W n -7 - /„_ 5 W„_ 6 

w n -i = M„_ 4 - 7„_ 4 w„_ 8 - a„_ 4 u;„_ 7 - .g„_ 4 w„_ 6 - f n -4W n -5 
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Cn-1 = d n -i - fci^i if c„_i = then c„_i = t. 
»=i 

Step6: i = 5, 6, . . . , n — 5 
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V n ~A = Kn-4 ~ 7n-4 w "-8 ~ an-4«n-7 ~ 9n-iV n -e — /n-4«n-5 

V n -3 = Af„_ 3 - 7 n _ 3 W„_ 7 - a n _3U„_6 ~ .9ri-3"n-5 - fn-3Vn-4 

Vn-2 = Al-2 ~ 7„-2 W "-6 ~ ttn-2"n-5 ~ 3n-2fn-4 ~ /n-2«n-3 
n-2 



ra-l 

c« = d rl — /ijWi if c„ = then c n = t. 

4=1 

n 



Step7: Compute det(X) = JJ Ci . If det(AT) = 0, then OUTPUT(Singular 

i=i 

Matrix) stop. 

Step8: For i — 1, 2, . . . , n compute and simplify the components Si jTl , Si tn —i, 
Si, n -2, Si <n -3, Si, n -4) S^ n -5 of the col umns Cj where j = n,n - l,n - 2,n - 
3,ra-4,n-5 by using (14I5I6I7I8I9I10I) . 

Step9: Compute the components of the columns C„_6 and C„„7 by using 
(|11I12[) . then for j = n — 8, n — 9, . . . , 1 and i = 1, 2, . . . , n compute and simplify 
the components S!y by using (TT3]) . 

SteplO: For i, j = 1, 2, . . . , n substitute the actual value t = in all Sij. 

Let R be an n x n matrix as in the following form: 



1 

1 

R 











It is clear that R is a nonsingular and its inverse matrix is itself. Let Y be a 
cyclic anti-nonadiagonal matrix. Since there is the following relation between 
the cyclic nonadiagonal and the cyclic anti-nonadiagonal matrices 

Y = KR, 

the inverse matrix of Y is obtained as 



Y- 1 = RK~ 
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3 Numerical Example 



Example 1 Consider the matrix D as in the following 
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We apply the algorithm to it and we have 



(ci,c 2 ,c 3 ,c 4 ) = (1,4,3, -|) 
(53,54) = (1,3) 

(ei,e 2 ,e 3 ,e 4 ) = (1,-1, -1,-^ 
{h,k 2 ,k 3 ,k 4 ) = 
(h 1 ,h 2 ,h 3 h 4 ) = (l,0,l§) 
{vi,v 2 ,v 3 ,v 4 ) = (1,3,-1,-1) 



(/2,/ 3 ,/ 4 ) = (-2,0,1) 
a 4 = 1 

(T U T 2 ,T 3 ,T 4 ) = (2,1,-1,1; 
(wi,w 2 , w 3 ,w 4 ) = (1,2, -1, - 
(P U P 2 ,P 3 ,P 4 ) = (-1,5,0,|) 
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(7 5 >7 6 ; 77>78;79,7io) = (-1,4, 
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det(lf) = 4715 
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In addition the inverse of the anti-nonadiagonal matrix Y is obtained as 
follows: 
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